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We discuss a recent approach for overcoming the poor convergence of the perturbative 
expansion for the thermodynamic potential of QCD. This approach is based on self- 
consistent approximations which allow for a gauge-invariant and manifestly ultraviolet- 
finite resummation of the essential physics of the hard thermal/dense loops. The results 
thus obtained are in good agreement with available lattice data down to temperatures of 
about twice the critical temperature. Calculations for a plasma with finite quark density 
(i.e., with a non-zero chemical potential //) are no more difficult than at \i — 0. 

1. The failure of the conventional perturbation theory 

The most compelling theoretical evidence for the existence of the quark-gluon plasma 
comes from lattice QCD which shows a clear signal for a deconfinement phase transition 
[[L|. Above it, the lattice results slowly approach the ideal-gas limit from below, with 
important deviations, though, of about 15-20%, up to temperatures ~ 5T C , and these are 
expected to remain noticeable (~ 10%) even at temperatures as high as 10 3 T C MM. 

This suggests a picture of the high-temperature phase of QCD where the interactions are 
more important than one would naively expect on the basis of the asymptotic freedom 
alone, but where the effects of these interactions are nevertheless small enough to be 
computable via weak-coupling techniques. 

The weak coupling expansion of the thermodynamic potential T (free energy) is present- 
ly known M to order a^ 2 , or g 5 (a s = g 2 fiir), but it shows a disappointingly poor conver- 
gence except for coupling constants a s ^ 0.05 (which would correspond to temperatures 
<; 10 5 T C ). Already the next-to-leading order correction of 0(g 3 ) signals the inadequacy of 
the conventional perturbation theory except for very small coupling, because, in contrast 
to the leading-order terms, it leads to a free energy in excess of the ideal-gas value. 

This blatant inadequacy is somewhat surprising, since one expects the free energy to be 
dominated by the hard thermal fluctuations with momenta k ~ T, for which perturbation 
theory should apply. Indeed, at temperatures T ^ T c ~ 300 MeV, the QCD coupling is 
reasonably small, a s ^ 0.3, when renormalized at the Matsubara scale ft = 2tcT. 

But a closer inspection of the perturbative expansion reveals that this is truly an ex- 
pansion in powers of g (rather than a s ), with g ~ 1 for all temperatures of interest, and 
that the largest "corrections" are associated with odd powers of g. The latter come from 
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resummations which take into account the phenomenon of Debye screening at the scale gT 
||. Thus, the large perturbative corrections are actually associated with soft degrees of 
freedom, with momenta of order gT, and arise when the contribution of the soft modes to 
the free energy is expanded in powers of g. But this expansion is potentially troublesome, 
since, as we shall shortly recall, the soft modes are non-perturbative. 

2. The quasiparticle picture of the quark-gluon plasma 

The soft degrees of freedom are collective excitations which would not even exist in 
the absence of interactions ||. To leading order in g, their dynamics is described by 
an effective theory obtained by integrating out the "hard" (k ~ T) thermal fluctuations 
to one-loop order. This generates a set of non-local self-energy and vertex amplitudes, 
known as "hard thermal loops" (HTL) M, which encompass screening effects and non- 
trivial dispersion relations. At momenta k ^ gT, the HTL's are leading-order effects 
and must be resummed for consistency ||. Thus, when computing the contribution of 
the soft modes to thermodynamical functions, these modes should be treated as dressed 
quasiparticles, with properties described by the HTL effective theory. This suggests a 
description of the thermodynamics of hot QCD in terms of weakly interacting "hard" 
and "soft" quasiparticles (rather than the more strongly interacting elementary quanta). 
This is also supported by the success of phenomenological fits involving simple massive 
"quasiparticles" in reproducing the lattice results 0. 

Quite generally, the physical information on the quasiparticles is contained in the spec- 
tral density p(u, k) related to the corresponding propagator by: 

d(w, k) = rppiR 

J-oc 2,71 k — UJ 

For free massless excitations, Po(uj, k) oc 5(u 2 — k 2 ). In the HTL approximation, the spec- 
tral densities are divided into a pole piece at time-like momenta, and a continuum piece 
at space-like momenta. At soft momenta k ^ gT, all the pieces of the HTL spectral func- 
tions appear to be equally important in our numerical calculations of thermodynamical 
quantities |§. At large momenta, the HTL spectral densities take the approximate form 

p(u, k) « 5{u 2 - k 2 - rr4) for k ~ T, (2) 

where m 2 ^ ~ g 2 T 2 is the leading-order thermal mass (or "asymptotic mass") of the 
hard excitations. Thus, quite remarkably, the HTL approximation describes correctly the 
mass-shell behaviour at both soft and hard momenta. 

In traditional perturbative calculations of the thermodynamics performed in imaginary 
time |J, the HTL's play almost no role: only the Debye mass ~ gT needs to be 
resummed in the static (uo = 0) electric gluon propagator [|J . Such a simple resummation 
retains only one moment of the spectral function in eq. (|I|). Although this is enough to 
obtain the leading order contribution, oc g 3 , of the soft modes to T , it is clear that this 
procedure mistreats most of the physical content of the HTL's. 

To overcome this limitation, two approaches have been recently proposed to perform full 
resummations of the HTL self-energies in the calculation of thermodynamical functions 
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In Refs. [§], this has been done by merely replacing the free propagators by the 
HTL-resummed ones in the expression of the free-energy of the ideal gas: 

F Q = l - Tr log D X — > Tutl = \ Tr log^ 1 + n OTL ) . (3) 

In principle, this is just the first step in a systematic procedure which consists in re- 
summing the HTL's by adding and subtracting them to the tree-level QCD Lagrangian. 
This would be the extension to QCD of the so-called "screened perturbation theory" [Q, 
a method which, for scalar field theories, has shown an improved convergence indeed, 
in two- and three- loop calculations. But in its one- loop approximation in eq. (|^), this 
method over-includes the leading-order interaction term oc g 2 (while correctly reproduc- 
ing the order-g 3 contribution), and gives rise to new, ultimately temperature-dependent, 
ultraviolet divergences and associated additional renormalization scheme dependences. 

Our approach on the other hand [§] is based on self-consistent approximations using the 
skeleton representation of the thermodynamic potential which takes care of overcounting 
problems automatically, without the need for thermal counterterms. 

3. The (approximately) self-consistent entropy 

Specifically, we consider the 2-loop self-consistent (or "^-derivable" [|(J) approximation 
to the thermodynamic potential J 7 , and focus on the entropy[|, which in this approximation 
takes a simple, effectively one-loop, expression [with N(ko) = l/(e /3k ° — 1)] : 

S = - J — — {imln/r 1 - lmU[D] ReD} (4) 

but in terms of fully dressed propagators so that U[D] is the one- loop self-energy built out 
of the propagator D. Thus, any explicit two-loop contribution to the entropy has been 
absorbed into the spectral properties of quasiparticles. The price to be paid is that Dyson's 
equation D^ 1 = Dq 1 + U[D] becomes an integral equation, which is further complicated 
by UV problems and, in gauge theories, also by the issue of the gauge symmetry. 

In spite of these complications, the expression (f|) has some obvious virtues: In addition 
to its simple quasiparticle interpretation (as the entropy of a non-interacting gas of quasi- 
particles with effective propagator D), it is manifestly ultraviolet finite (the derivative 
of the statistical factor acting as an UV cut-off), and provides a non-perturbative ap- 
proximation to the thermodynamics which is perturbatively correct up to, and including, 
0(g 3 ) (since the neglected 3-loop diagrams start contributing at 0(g A )). 

To cope with the problem that ^-derivable approximations are not gauge invariant 
in general, we have proposed gauge-independent but only approximately self-consistent 
dressed propagators as obtained from (HTL) perturbation theory. Using these in eq. (|]) 
gives a gauge-independent and UV-finite approximation for the entropy, which, while 
being nonperturbative in the coupling, contains the correct leading-order and next-to- 
leading order effects of the interactions. Both arise from kinematical regimes where the 
HTL's are justifiable approximations — at hard momenta, they involve the HTL's only in 
the vicinity of the mass-shell, where the HTL approximation is sound, cf. eq. (0). 

^^More generally, on the first derivatives of the thermodynamic potential, like the entropy and — for 
plasmas with non-zero chemical potential — also the quark density. 
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Figure 1. Entropy and pressure for pure-glue SU(3) Yang-Mills theory: Approximately 
self-consistent resummation [Q vs. dimensional reduction || and 4-d lattice data 



In the left half of Fig. |T], the results for the entropy in a pure HTL approximation 
and in a next-to-leading one (NLA) which includes soft corrections to the hard asymptotic 
mass are compared with lattice data for pure-glue SU(3) theory, showing good agreement 
for T ^ 2.5T C ; in the right half, the resulting pressure?] is compared in addition to the 
calculation of the pressure in dimensional reduction of Ref. || for temperatures up to 
10 3 Aq CD . The fairly good agreement with the best estimates of Ref. || (e = 10) seems 
to validate our assumption of comparatively weak residual quasiparticle interactions. 

Our method has been applied successfully also to plasmas with non-vanishing quark 
density (i.e., non-zero chemical potential) [Qj, for which lattice results are not yet available. 
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2 with renormalization scheme adjusted so as to match that of Ref. 



